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This work is intended to be a contribution to the study of the morphology of the 
rising convective columns, for a better representation of the processes of entrainment 
and detrainment. We examine technical methods for the description of the interface 
of expanding clouds and reveal the role of fingering instability which increases the 
effective length of the periphery of the cloud. Assuming Laplacian growth we give a 
detailed derivation of the time-dependent conformal transformation that solves the 
equation of the fingering instability. For the phase of slower expansion, the evolution 
of complex poles with a dynamics largely controlled by the Hilbert operator (acting 
on the function that represents the interface position) leads to cusp singularities but 
smooths out the smaller scale perturbations. 

We review the arguments that the rising column cannot preserve its integrity (seen 
as compacity in any horizontal section), because of the penetrative downdrafts or the 
incomplete repulsion of the static environmental air through momentum transfer. 
Then we propose an analytical framework which is adequate for competition of two 
distinct phases of the same system. 

The methods exmined here are formulated in a general framework and can be easily 
adapted to particular cases of atmospheric convection. 
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I. INTRODUCTION 


In the atmospheric convection leading to cumulus clouds there is a continuous exchange 
of heat and water (vapor, liquid) between the rising column and the environmental air. This 
problem is complex and has been examined by observation, analytical theory and numerical 
simulatioiP,^,!^,^,^. The diversity of situations cannot be captured by a unique model. 

The rate of exchange depends, besides many other factors, on the geometry of the contact 
interface between the two gaseous media. Part of the interaction with the environment takes 
place at the periphery of the cloud, an interface that evolves and is subject to geometric 
instabilities. In the expanding phase the small scale structure is determined by instabilities 
of the fingering type. Like in the case of the Hele-Shaw instability, one can map the physical 
problem on a time-dependent conformal transformation to a hxed complex line. At small 
scale, there can be random fluctuations due to the background turbulence. The effective 
length of contact between the cloud and the environment is then much greater than the 
perimeter of hull of the cloud in horizontal plane. This enhances the transport between the 
two neighboring gaseous media. 

In the phase where the expansion slows down (the input from the rising flux is reduced) 
the small scale prohles (resulted from “hngering” and random fluctuations) are smoothed 
out by a process of attraction and alignement of the complex poles that dehne the solution. 
The result is that the interface exhibits singularities of a special type, cusps. They occur 
through coalescence of small scale quasi-singular perturbations of the wrinkled interface. If 
the time evolution leads to self-intersection of the interface, a parcel of environmental air is 
simply swallowed into the expanding cloud. 

Finally, inside the rising column there is interaction between the cloud air and the en¬ 
vironmental air that either remains by breaking the initially compact rising column, or 
penetrates from above as downdrafts on long vertical distances. 

This last case requires few comments. There are various forms of contact between the 
two gaseous media that are in relative motion. In a simplihed descriptive picture we can 
see the rising column as a stream of gas with specihc properties penetrating a volume of 
stationary gas (static environmental air) with different properties. If the speed of the stream 
is high, there is transfer of momentum from the stream to the the environmental air through 
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collisions at the front of contact, i.e. at the top of the rising column. Then the environmental 
air is deflected (pushed up and sidewise) and the body of the column preserve its properties. 

Alternatively, if the linear momentum of the rising column is small, then the stream of 
the convective air cannot repulse through momentum transfer the environmental air and 
there is easier interpenetration of the two gases. The rising stream is teared apart and 
elements of the air of the column are interspersed between elements of the environmental 
air and this occurs up to small spatial scales, of the turbulence. Obviously, the exchange of 
heat and water (vapor and liquid) is much more efficient in this case. This is the case of 
either shallow convective events that are dissipated before becoming a buoyant column or 
(and this is most interesting) of clouds in the last phase of their ascent where the vertical 
advancement is slowed down and elements of the cloud are dispersed in the environment. 

For a comparison, the two cases have been found experimentally in the expansion of a 
laser-generated plasma in a low density plasmsP. With variation within a range of parameters 
the interaction changes from interpenetration to formation of an interface. 

These two situations are limiting cases. In general it is expected that the rising convection 
columns that form the cumuli are characterized by strong buoyancy that allows them to 
reach high altitude. Therefore they are closer to the first case described above, where the 
stream is able to push off (out of) its way the environmental air and the contact with this 
one mainly takes place at the moving top and at the circumference. Even if the vertical 
momentum of the rising column is high, the possibility to retain a compact structure is 
improbable. More realistic is the expectation that the initially compact column will break 
up into streams separated by irregular but vertically-connected volumes (channels) of static 
environmental air. In short, the breakup of the rising column leaves open spaces inside the 
initially compact area and in these spaces there still is static environmental air. The transfer 
of mechanical momentum from the rising cloud air to the static environment is mediated 
by exchange of turbulent eddies, therefore implicitely involves dilution of cloudy air and 
evaporative cooling for the environmental air. This occurs mainly at the top of the rising 
column. Parcels of the cooled environmental air will descend and will push the still un-mixed 
environmental air present in the contiguous channels remaining inside the broken column. 
The mass conservation requires the environment to respond to the rise of the convective air 
by currents of descending air, which takes place in a layer around the rising column but also 
in the spaces left open after the breakup of the rising column. Therefore, elements of rising 
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air that are at a certain level height will be in close contact with environmental air that, 
actually, originates at levels of height that are much above the current level and are not yet 
mixed with cloud air. Then the mixing will take place at various heights. This picture is 
very close to what is observed in a classical Reyleigh-Benard system in the transitory phase 
where the purely conductive state is going to be replaced by the convective state. This 
bifurcation is preceded by emission of thermal streams (plumes) at the hot plate, which do 
not have the chance to produce a full scale convection. However they are able to determine 
thorugh mass and momentum conservation, generation of opposite streams, originating at 
the cold plate and descending. There is no mixing between the rising and descending plumes 
except at late stages, when the convection span the whole volume. Parcels of environmental 
undiluted air, of large dimensions (up to 500 (m)) have been observed at all levels in the 
cloud. 

This is the situation that we have in mind when we consider the possible downdraft of 
environmental air not yet mixed with cloud air. A parcel of environmental air can be pushed 
to descent by the request of mass and momentum conservation, to compensate the rise of 
streams of rising cloud air, in a structure of the broken column characterized by coexistance 
of buoyancy-driven ascent and vertical channels of environmental air. We therefore note 
that the latter may be not yet cooled by the evaporation of the water after mixing with 
cloud air close to the top. 

Therefore we propose to include in the physical picture the downdraft of pure environ¬ 
mental air, not yet mixed with the cloud air, resulting from the constraints of mass and 
momentum conservation. The downdrafts are located inside the broken cloud column. In 
addition, after mixing, the new downdraft will be cooled by evaporation and the mixed 
parcel will descent even more, up to the cloud base. 

During the rise there is a smooth change from the situation of strong stream, specihc to 
the initial stages and the soft slowing down, characteristic of the phase where the column 
reaches the highest level. In the first phase the contact between the convective column and 
the environment takes place at the top, at the interface with non-mixed downdraft and 
at the periphery and in the last phase there is less momentum transfer and easier mutual 
interpenetration of parcels of air. The exchange of heat and water is less efficient in the first 
phase and is highly efficient in the last phase, which actually accelerates the process of loss 
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of buoyancy. 


This approximative and descriptive picture suggests to treat separately the two situa¬ 
tions. In the first one, the exchanges between the convective column and the environment 
require to examine the periphery. Close to the hnal slowing down of the ascending column, 
the exchanges involve a large volume where parcels of similar dimension of cloud and envi¬ 
ronmental air are intermingled. Since the top of the column has been in contact with the 
environment all along the rise, one must respresent the breakup of the rising column into 
smaller columns and the presence, between them, of columns belonging to the initially static 
environment or downdrafts resulting from evaporative cooling of parcels of mixed air. 

These two situations will be our subject. For the problem of cloud periphery, we will dis¬ 
cuss the possible evolutions of the interface between the expanding gas and the surrounding 
air: the shape of the interface, generated the hngering instability and random fluctuations. 
As expanding front, we examine the formation of cusp singularities. 

For the problem of breaking of the rising column, we will discuss the statistics of phase 
competition. We hnd that, in this schematic representation, based essentially on geometric 
aspects of competition of two distinct phases of the same gas, it is relevant to discuss in 
terms of labyrinth structnred in the horizontal plane. 

Therefore we have simultaneously a problem of separation of phases and a problem of 
interface dynamics. For tractability we divide the problem into two different components: 
the interface dynamics for each component, as well as for the large bonndary circumscribing 
the full horizontal section of the convective column, is investigated with methods of pole 
dynamics and/or wrinkled advancing fronts; the phase separation is treated with method of 
coupled lattice maps, in order to follow easily the breaking of the compact raising column; 
and, separately. 
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II. GEOMETRY OF THE INTERFACE BETWEEN THE RISING 
COLUMN AND THE ENVIRONMENT 


A. The exterior interface 

1. Introduction 

It is considered that the the circumference of the cloud at a hxed altitude ^ has an im¬ 
portant role in the exchange of heat and water vapors with the environment. The turbulent 
diffusion sustained by random eddy exchanges introduces environmental air into the mass 
of the cloud and there the exchange of heat / vapor modifies the tendency of the column to 
rise. During the phase of rise the edge of the convective column has a vertical motion relative 
to the static environment, which creates a layer of vorticity at the interface. There are two 
mechanisms that can affect this interface, depending on the relative velocity and viscosity. 
The Kelvin-Hlemholtz instability can get a positive growth and peripheric elements of rising 
column are rolled up to create the known “cats-eye” pattern in a vertical-plane section, i. e. 
a ring vortex with vertical principal axis. Alternatively, a combined lateral expansion and 
rise of the column can produce the ring vortex at the head of the column, i.e. due to the 
inertial resistance of the static environment, a typical moshroom head (well known from 
Rayleigh - Taylor instability). In both cases parcels of environmental air are absorbed and 
entrained being surrounded by cloudy air, thus facilitating the mixin^^^. 

If at the circumference the transport processes associated with the entrainment depend 
on the area of contact between the cloud and the environment then a careful representation 
of this area is necessary. Equivalently, at hxed altitude z, a good representation of the 
geometry of the interface cloud/environment is necessary. We will discuss small scale and 
respectively large scale structure. The small scale structure of the interface is generated 
by two processes: (1) deterministic instabilities, like fingering; (2) random perturbation, 
related to the turbulence. On a large scale, the structure of the interface can lead the cloud 
to incorporate (swallow) parcels of environmental air. The mechanism originates in - and 
is a limiting form of, - the cusp singularity that is formed as an asymptotic organization of 
the small scale quasi-discontinuities. 
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2. The small scale structure: the “fingering” instability of the interface 

The interface between the air of a rising column and the environment air shows a specihc 
prohle. In every horizontal plane there is a fluid (cloud air) expanding into a static fluid 
(the environment). An universal model for such process assumes that the velocity of the 
interface is derived from the gardient of a scalar function that verifies the Laplace equation 
(is harmonic function). Then the interface is subject to ^fingering” instability. The role of 
the Laplacian field is played by the pressure of the expanding gas. We need an analytical 
instrument to describe the small scale breaking of the continuity of the derivative of the 
line of the interface. This wiggled prohle is the place where the exchange of heat and water 
vapor takes place. 

The two huids (cloud and external environment) are assimilated with two different 
phases separated by a moving interface T (t), a curve in the physical plane of coordinates 
[X [t) ,Y (t)]. In a simplihed representation, the interface is a line that extends between 
—cx) and +cx), i.e. it separates two regions of the plane. The region I is the inside of the 
cloud, limited by T [t) and the region II is outside, the environment. It is assumed that the 
velocity of the expanding huid (the cloud) is the gradient of a scalar function P (X, Y). 

v„ = — (VP)^ at the interface T (1) 

The subscript n means projection of the vector VP on the normal at the interface T. The 
equation for the scalar function is 2D Laplace: AP = 0 in the cloud region, he.in the lower 
part limited by T (t). This function is P = 0 in the free (environment) region 11. The physical 
source of expansion is the input of cloud air from below the current position of the interface 
r(t). This is represented as an asymptotic condition for velocity: somewhere very far inside 
the cloud (Y —)■ —oo), the velocity of the air is a constant directed toward the interface 

VP = ey for V —)■ —oo (2) 

In addition 

P = 0 at the interface T (3) 
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3. The time-dependent complex conformal transformation 


The evolution consists of changes in time of the curve T (t) = [X (t ), Y (t)] representing 
the interface, i.e. expansion of the boundary cloud/environment in any horizontal section of 
the cloud. The idea is to hnd a mapping between the physical plane {X, Y) and the complex 
plane z = x -\- iy. At every moment of time t, the lower semi-plane in the mathematical 
complex plane {x,y) is mapped to the space below the interface T (f), where P verihes the 
Laplace equation. The evolution of the interface is then a set of conformal transformations 
parametrized by timd^. The scalar function P, is defined as the real part of a new complex 
variable, whose imaginary part is a function {Z) 

W{Z) = P (Z) + (Z) (4) 

The two real functions T and P are harmonically conjugated where W is holomorphic. 

The conformal map is the function /, / : C — )■ C, Z = X + iY = / {z, t) where 2 ; = x+iy. 
Since / (z, t) maps the lower half complex plane {y < 0) to the region under T (t),Y <0, its 
derivative ^ should have no singularities or zeros in the lower semi-plane. All of them must 
be in the upper semi-plane. Translating Eq.(|^ it is expected that at very large distances 
on the plane, relative to the interface, the variables will have very close values 

fiz,t) 2 : for 2 : —)■ X — ioo (5) 

which corresponds to constant velocity of the incoming air, the source being the air rising 
from below the horizontal plane of the current height. 


Now the system is rewritten: the new function (instead of the pressure) is W {z) and the 
new variables (instead of (x, y) ) are 2 : and 


dW 

dz 

dW 

dz 

ReW 


0 {W is holomorphic) 

—i for .2 —)■ —ioo at large distance from T 
0 at 2 ; = X — iO (at the mapped interface) 


( 6 ) 


From Eq.([^ 


Re 


n 
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0 


(7) 









here n is a complex number associated to the normal versor at the interface. The solution is 

W (z) = -iz (8) 

from where the scalar function P is 


F(X,F)=Im [f-^{X + iY,t)\ 


The normal at the interface is 


n = —i 


dz 
dz \ 


at z = X — iQ 


Then the equation describing the Laplacian growth is 

= 1 at z = X — iO 

(for z on the real x axis, just below) 


Im 


df {z,t) df {z,t) 


dz 


dt 


(9) 


( 10 ) 


( 11 ) 


which is the “Polubarinova - Galin” equatioiP,®,™™. 

A particular form of the solution / is derived by Ponce Dowson and Mineev Weinstein in 
Refspra 

Since / is holomorphic in the lower complex half-plane a general expression is 
dehned by choosing for ^ a number of zeros and poles in the upper half-plane. This provides 
the explicit form of the mapping, at a hxed moment of time. Since we have assumed inhnite 

extension of the interface, a class of solutions is 

v+i 

/inf {z,t) = z - it log [z - (i (t)] (12) 


1=1 


where 


and {a', a”) are real. 


ai = a[ + ia'l are iV + 1 complex constants 


Q = + irji are iV -f 1 singularities 


( 13 ) 

( 14 ) 


simple poles of ^ that move in time. 


A detailed treatment of the conformal mapping is provided in the Appendix A. The 
evolution keeps the zeros and the poles in the upper half-plane, where they are initialized. 


At any moment of time the new positions of the singularities are inserted in Eq.(12) and 


the conformal mapping is determined. Therefore we have the new shape of the interface 
between the two physical media (expanding cloud and environment). 
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Time evolution of few singularities ^(t) 



Figure 1: The curves represent the time evolution of the positions of the iV + 1 singularities 
Ckit) = [^kit), rjkit)], {k = 1, N + 1 for N = 9). The initial positions rjkit = 0) are at in a 
small interval around 10 and the real parts = 0) are distributed on equal intervals. 
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X=Real(f) as function of XL 



Figure 2: The lines X (t) = Re [/ {z,t)] at 5 moments of time, t = 90,..., 95. Since both 
real and imaginary parts of Q are not too close to 0, the lines are smooth. 


Y=lmag(f) as function of XL 



Figure 3: Same as Fi^for Y (t) = Im [/ {z,t)]. 
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The shape of the Interface: Y=lmag(f) as function of X=Real(f) 



Figure 4: The interface at 5 moments of time, as in Fig^and]^ 
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B. The small scale structure: roughening of the interface due to random 
fluctuations 


Besides the deterministic evolution described by the fingering instability, one should 
consider random fluctuations that produces the roughening of the interface. It has been 
found that the interface roughness in 2D is algebraic on short length scale. Here we only 
mention from Ref.™ a part of the argument that introduces the Hilbert transform. This 
is an analytical step which prepares the discussion on the large scale features of the cloud 
boundary. We consider the interface between the cloud and the environment consisting, in 
the horizontal plane, of a line T of length L. The coordinate along the interface is s and 
along the local normal is y. The position in plane of the current point on the interface is 
given by the distance 

h (s, t) 

relative to a fixed reference system. The Laplacian held (the pressure) cj) (x, t) is defined as 

= 0 for y < h inside the cloud, and (15) 

0 = 0 for y > h 


The velocity of a current point of the interface is given by the gradient of the scalar function 
0, as 


dt \dy ds ds) 


(16) 


This equation describes the gradient flow. The last term describes the local dilation or 
the compression of the length of the curve T. The deterministic part of the dynamics is 
introduced by the average motion of the front 


h = Vt 


A small perturbation of the “height” h {x, t) taken as 


hq (s) = hq (t) exp (iqs) 


acting on the moving interface decays as 


ft, (s) ~ exp |-ff ((/) (] 


( 17 ) 


(18) 


(19) 
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so the perturbation is exponentially vanishing with the rate 


a{q) = V\q\ (20) 

Then the relaxation of the interface “height” is a simple linear decay of the logarithm of hq, 
working in the Fourier space. Since we want to take into account the random fluctuations, 
it is introduced a noise source 

=-V\q\hq{t)+r]q{t) (21) 

This is a nonlocal equation since a {q) acts in Fourier space. The noise is Gaussian 

{Vq (t) Vq' {t')) = jH(l + Q') t') (22) 

We note that the dynamical equation for the noise-driven interface contains the term 
—V\q\hq{t). It comes from taking the Fourier transform of the real-space function h, 
multiplying with the absolute value of the Fourier space variable |g| and returning to real 
space. This is the Hilbert transformation applied on h and will play an essential role in the 
following. 

In Ref.l^ it is shown that the width of the interface increases as the log of the length L. 


C. The large scale dynamical structuring of the interface: cusp singularities 


We are now interested in the phase of the cloud expansion where the convection flux 
coming from lower levels is progressively reduced, the column reaching a regime of quasi- 
stationarity. On a large spatial scale (of the circumference of the expanding cloud) the 
interface has a dynamics of the expansion of a front as the propagation of a planar flame 
into a static homogeneous medium. A model equation for the latter case has been developed 
by Sivashinsk}^,^. It treats the advancement of a front of a flame in a chanel. 

The planar flames expanding freely has an interaface that is unstable. For a simple model 
in ID (flames propagating in a channel of width L) the variable is 


h [x, t) = position of the flame front above the x axis 

The equation of Sivashinskj^ is 
dh{x,t) 1 

Ft 2 


dh {x, t) 
dx 


. r, / M d^h (x, t) 
A{h{x,t)} = u +1 


(23) 
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in the domain 


0 < X < L 

The functional A is the Hilbert transform. To dehne its action, hrst one makes the Fourier 
transform of the function h {x, t), 


h (x, t)= exp (ikx) h {k, t) dk 
J —OO 

then multiply h {k, t) by the absolute value of the Fourier variable 

A{h{k,t)} = \k\h{k,t) 


( 24 ) 


( 25 ) 


and returns to the real space, A {h (x, f)}. We can trace the occurence of the Hilbert operator 


in the equation for the front advancement from the derivation of Eq.(21). 


The front is unstable and it generates singularities in hnite time. The singularities are 
of the cusp type. We note that this evolution implicitely renders the interface piecewise 
smoother since it collects the smaller scale singularities into a single, giant cusp. This 
corresponds to the late phases, where the flux of cloud air is reduced and the reserve of 
buoyancy is decaying. 


1. Formalism: Pole decomposition 

The origin of the cusp prohle can be understood from the evolution of the positions 
of the complex singularities (poles) of the solution. The motion of the pole singularities is 
controlled by the Hilbert operator. The multiplication with the absolute value of the Fourier 
variable \k\ is equivalent to the differential operator It makes the poles to approach 
the real axis, from both sides. The accumulation of poles (“condensation” at the limit of 
continuous density of poles) leads to the formation of the cusp, a singularity of the interface 
that is qualitatively similar with what is frequently seen in the evolution of an expanding 
cloud. This suggests that the large scale structure at the expansion of the circumference of 


a cloud as it rises can be studied using the Sivashinsky equation, (23). 


Therefore the hrst step in an analytic development is to extend the space coordinates to 
complex variable^,^,^,^,^,^. The nature of processes that take place during time evolution 
is: (1) attraction between the poles along the horizontal axis (x or 6) leading to clusterization 


15 




of poles along vertical direction in the complex plane, and (2) dynamical evolution of the 
poles towards the real axis. This produces singularities that look like cusps. 

The problem is restricted to a single space variable, relative to which we measure the 
spatial position h {x, t) of the interface. The nature and the positions of the wrinkles of the 
interface can be associated to the singularities of the interface - function in the plane of the 
complexihed spatial variable. A singular prohle will occur in hnite time but for short time 
the equation is integrable which means that it has the Painleve property and the formal 
solution is meromorphic and can be written as an expansion in a set of order-one poles. 

In Ref.^it is porposed the Lagrangian version of Eq.(21) with the white noise converted 


into a diffusion generated by viscosity u, essentially constructed on the ground of Burgers 
equation 

rl^i Flu Flu'll 

( 26 ) 


du du ^ . d‘^u 


As before, A [m] : u{k,t) —)■ \k\u{k,t) where the Fourier transform of the function u has 
been introduced M (x, t) = dk exp{ikx) u{k,t). The solution is a meromorphic function 
expressed in terms of simple poles by terms like 


Pa = 


X 


(27) 


The operator A applied on such term results in 

Apa = A ' 


X 


(28) 


d 


= sign [Im (za)] i-^Pa (x) 


The operator A produces advection of the poles, along the imaginary direction, towards the 
real axis (both from above and from below since the operator multiplies with |fc|). This 
enhances the effect of the singularity, i.e. the real function which is the solution becomes 
even more perturbed there. The typical prohles are quasicusps, which are cusps with the 
tip rounded due to the presence of u. In Refs.^,^ this is explained by showing that the 
singularities never reach the real axis, but remains at a distance whose expression is given 


by z/. This will produce wrinkles. The solution of the Eq.(26) can be written, for the linear 
geometry as 


2N 


U 


{t,x) = -2z/^ 


^ X - (t) 


(29) 
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where Za (t) are 2N poles placed symmetrically, as complex conjugated pairs {u (x, t) results 
real) relative to the real axis. The equations of motion is 


dza 

dt 


I3^a 




isign [Im {za)] 


( 30 ) 


The poles tend to attract themselves horizontally, parallel with the real axis. As noted 
in Ref.l^ this results from the dominant behavior of the x coordinate of the poles 


dxj 

dt 


-J^sin {Xj-Xk) [...] 


In the paranthesis [...] there are terms that are invariant to the change xj O Xk, for any pair 
j, k = 1,N, with j 7 ^ k. This factor is positive. When 0 < Xj < Xk we have dxj/dt > 0 and 
for 0 < Xfc < Xj we have dxj/dt < 0. In both cases, Xj evolves to be closer to Xk- 

In addition there is a tendency of the poles to place themselves on a line parallel to the 
imaginary axis (this is the analog of the same phenomenon for the Burgers equation). At 
short range the poles tend to repel each other vertically and the repulsion between two poles 
aligned on a vertical {y) line becomes inhnite when they come too close. At longer distances, 
the interaction become attractive. 


Details on the integration of Eq.(30) are in Appendix B. 


2. Further developments 

The stability of the interface consisting of a giant cusp has been examined irP^ ancP^ 
where it is explained the role of the noise in generation of new poles in the structure of 
the solution. This may also explain the self-fractalization of the interface resulting from 
self-acceleratioiP of the flame front determined experimentallj^^. 


III. THE BREAKING OF THE RISING CONVECTION COLUMN 

A. Introduction. The loss of compacity of the column 

During its rise the convective column may undergo breakup and loss of compacity in its 
volume. The column will then consist of vertical streams of cloudy air separated by volumes 
of environmental air, also vertical and of irregular shape, which is either static and not yet 
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Time evolution of the singularities ^(t) 



Figure 5: The time evolution of the complex poles of the solution of the interface problem 

for the case of the cusp singularity. 

mixed or consists of downdrafts after evaporative cooling of parcels entrained close to the 
top of the column. As recalled in the Introduction, the simplest view on the convection is a 
picture of intepenetration of two gases, one being the stream of rising air whose ascension 
is sustained by buoyancy and the other being the environment. There is a wide spectrum 
of possible evolutions, from the early loss of the streaming compacity to the preservation of 
robust columnar features during all rise. We try to hnd an analytical representation of the 
breaking of the rising column. 

There are observations revealing that, on a horizontal plane, the updraft inside a cloud 
presents substantial inhomogeneit}^,l^,l^. There are large variations of the vertical velocity, 
in magnitude and sign, meaning that there are updrafts and downdrafts and also regions 
of static environmental air. It results that there is a breaking-up of the rising column. 
Obviously, the surface of contact between the convective air and the environment, either 
static or penetrative downdrafts, is much larger than in the case where the column remained 
compact as it rises. The examination of the breaking up of the column, with substantial 
increase of the “interface”, is a necessary step in a better representation of the exchanges 
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The profile of the interface function h(x,t) 
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Figure 6: The profile of the interface at several moments of time. The initial state has 
many oscillations of small amplitude and the hnal state presents the cusp singularity. 


and transport processes. In particular this refers to the model of Squired^ where it is 
assumed that the environmental air is mixed with the cloud at the cloud top. Then, due 
to evaporative cooling, the mixed air loses its buoyancy and descends deep into the cloud 
column (kilometers 1^. The exchange of heat and vapors with the cloud air continues for 
these internal downdrafts and the result is dilution of the cloud air but also increase of 
the buoyancy of the air penetrating from above. The equilibrium between the penetrating 
downdrafts and the cloud air surrounding it is reached at equal buoyancy. 

We are interested in the process of competition between the rising air of the convective 
column and the environmental air. The result is the physical breakup of the column and the 
coexistence, at every horizontal level, of regions of rising air and of environmental air. The 
inhomogeneity of the cloud can be pronounced: inside the clouds there are narrow updraft 
regions but between them there are strong negative vertical velocity flows, due to rapid 
downdraft^. This pattern, in the horizontal plane, justifies the use of a model of interface 
dynamics which exhibits the labyrinth instability. 


19 




8 

6 

4 

2 

0 
x" 

13 

-2 

-4 

-6 

-8 

-10 

0 0.2 0.4 0.6 0.8 1 1.2 1.4 

xl 

Figure 7: The slope of the interface. 

B. The breaking of the convective column as a phase-competition dynamics 

We look for a schematic analytical and/or numerical description to the real phenomenon 
of breaking of the rising column into distinct vertical streams separated by regions of environ¬ 
mental air. It is easier to restrict to horizontal planes. The most elementary representation 
consists of the separation of phases of a fluid in 2D. For a binary fluid the variable is the 
concentration c where the two pure phases have c = ±1. The dynamics is described by a 
parabolic equation where the local change of the concentration is the Laplacian flow of the 
density of a functional of c, which can be called “chemical potential”. The free energy decays 
to zero when there is full separation of phases, i.e. the two phases occupy disjoint regions 
in plane. These regions can look similar to a labyrinth patterrP^^. For example, in the hnal 
stage of the slowing down of the rise of the cloud, where substantial loss of buoyancy has 
resulted from mixing with environment, the horizontal plane is mostly covered by regions of 
environmental air, with few patches of cloud. 

Technically, we can use the analytical approach based on one of the standard models, 
in particular Cahn-Hilliard system, studied in Ref.^, or cellular automataP. However it is 
more useful to implement another aspect besides the phase-separation: phase competition. 


The profile u{x,t) = dh/de 
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Time evolution of the profile of the interface 



Figure 8: The relief view of the evolution of the interface. 


In similar cases it has been adopted a representation of dynamical phase transition which 
employs a system of coupled cubic map lattic^. 

If the two phases have the same stability then the motion of the interface is governed as 
in a gradient flow, by the local curvature. If the two phases have distinct stability properties 
then the most stable phase will advance irreversibly into the other and will eventually replace 
it. The problem belongs to the same class as phase separation and nucleation phenomena. 

The coupled lattice maps are described by the set of equations for a discretized variable 
q (x, t) which represents the nonconserving order parameter at the points of a regular lattice 




flQiht)] 


+7 


j=n.n 




(31) 


where i = •••) is a set of integers that specifies the position of a point in the lattice. 

In our case the dimension of the problem is d = 2 (plane), i = (ii, * 2 ) = {ix, iy)- The sum in 
the square paranthesis extends over the nearest neighbors [n.n) of the point i of the lattice. 
They are in number of p in general and in d = 2 they are p = 4. The square paranthesis 
is actually a discretization of the Laplace operator and this term represents the diffusion. 


21 
















The nonlinearity of the dynamics of the order parameter q (x, t) is introduced by the term 

f[(l{ht)] = -q^ + {l + e)q + c (32) 

We recognize easily the meaning of choosing this nonlinearity. The “potential” / has two 
extrema, corresponding to non-dynamical equilibria. The order parameter can take one 
or another of these two equilibria values, and they are associated to the two states. The 
dynamical equation for q (x, t) is a discretized form of the Landau-Ginzburg equation. 



Figure 9: The end stage of the competition between two phases (which we would attribute 
to cloud, respectively environment). The vertical axis is the height of successive planes 


that are labelled by the time in the iterative map Eq.(31). The parameters are a = 1, 

e = 0.65, c = 0.10 and 7 = 0.09. 


C. Phase separation for equally stable states 

This corresponds to 

c = 0 


(33) 
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Figure 10: The same as Figj^ for a = 1, e = 0.45, c = 0.6 and 7 = 0.06. 


and the initial state is 


qo{i,j) = 0 (stable) (34) 

(zgomot) 

The spatial competition produces clusters of one phase inside the other. Any curved 
frontier between the two phases evolves with a velocit}^ 


V = —jK (35) 

transversal on the boundary. K is the curvature of the boundary. For a disk of radius R 
the curvature K = R~^ and the evolution is symmetrical with 

R^ (t) - R^ (0) = -27t (36) 

The initial situation is a compact covering of the domain of interest with a single phase. 
This initial state is perturbed by the random nucleation of the other phase on small regions. 
The evolution consists of extension of the regions associated with the Phase II, against the 
Phase I and the break up of the region initially occupied by only the Phase I. The process 
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0 100 

Figure 11: The competition of the two phases is shown as results from iteration. The 
structure (barely visible) has the small scale oscillations related to the propagation mode 

at the level of one mesh cell length. 


leads to the loss of compacity of the region of Phase I (cloud) and to progressive reduction 
of the surface occupied by it in horizontal planes. 

It has been provecP^that the iterative update adopted in Eq.(31) and (32) is equivalent 


with the curvature flow Eq.(35). The dynamic structure function, defined as the discrete 
Fourier transform of the correlation of the held q has at large time a decay like 

We note however that this is a purely geometric property, and the cloud-environment mixing 
is not taken into account. We have however a lower bound to the rate of disappearence of 
compact parcels of cloud in the late phases of the convection. 
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Figure 12: Left: Same as Figj^ after artificially removing the small scale oscillation. This 


makes more visible the two phases. Right: The contour slices of the left hgure, that are 


used for calculation of the length of the circumference of horizontal patches. 
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(g) 

Figure 13: The contour levels corresponding to patches of compact cloudy air in several 

horizontal planes. 













(a) (b) 

Figure 14: Area and length of the perimeter of compact patches of cloud air, with 
dependence on the height of the horizontal plane. The latter variable is mapped from the 
variable time of the iterative sequence. The parameters are: a = l, = 0.45, c = 0.6, 

7 = 0.06 and time step 0.02. The hgures correspond to 200 time steps. 
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IV. CONCLUSIONS 


We have examined three geometrical aspects that can be important in the quantitative 
studies of the exchanges of heat and vapor between the cloud and the environment. 

The first model, regarding the fingering instability can be an important step in repre¬ 
senting the fractalization of the boundary of the clouds. 

The second model, intended to allow a quantitative description of the cusp singular¬ 
ity of the interface, can be also implemented in any study of the lateral exchanges cloud- 
environment. It can be developed further for the study of the process of absorbtion of 
parcels of environmental air inside the cloud, but this is a problem of complex function with 
a certain difficulty. 

Finally we have considered the possibility to describe the loss of compacity of the rising 
column by a discrete, coupled lattice map model. At least the basic facts of decay of the 
convection and loss of continuity of the cloud tower, can be examined using this model. 

A wide range of similar models can now be proposed and further development can be 
considered. 
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A. APPENDIX. THE INTERFACE DEVELOPING FINGERS IN 
LAPLACIAN GROWTH 

A. Interface between infinite regions 


The evolution of the front of the expanding cloud, in any horizontal plane, is repre¬ 
sented through the conformal mapping from the lower half of the “mathematical plane” on 
the region below T, the interior of the cloud. The Laplacian growth is described by the 

Polubarinova - Galin (P-G) equation 


Im 


df {z,t) df 

dz dt j 


1 at z = X — it) 


(A.l) 


(for just below the real x axis) 



A class of solutions 


N+l 

/inf {z,t) = z - it - I'^ai log [z - Cl (f)] 

i=i 


(A.2) 


where 


and (a', a”) are real. 


ai = a[ + ia'l are iV + 1 complex constants 


Ci = Ci + i^i cire iV + 1 singularities 


simple poles of that move in time. 


(A.3) 


(A.4) 


Let us calculate explicitly 


dl 

dz 


N+l ^ 

1 —i > ai - 


and 


dl 

dt 


1 dCi it) 
^ - (lit) dt 


Ai+l 


df . . _ 1 


dCi (t) 


z-Ci (t) dt 


and the product is 

df^ 

dz dt 


Af+l 


Ai+1 


dCi D 




= I 


N+l 

1=1 

N+l 


z-Ci (t) 

1 dCi (t) 


+iy ai- - 

tr dt 

Ai+l A^+l 


(A.5) 


(A.6) 


(A.7) 


(A.8) 


1 _ 1 dCk jt) 

This must be taken for the line that in the mathematical (complex) plane z represent the 
interface 

z ^ X — iO (A.9) 
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which means that z will be replaced everywhere with real x. Before doing this we make a 
test, looking at what this expression becomes when we replace 


t (^rn (t) 


(A.IO) 


and find for the first two terms 

N+l 


N+l 


I ^ I ■ it) 

N+l ^ N+l ^ 


(A.ll) 






dCi (t) 


tr Cm{t)-Ci{t) tt Cn.{t)-Ci{t) dt 
According tcP,^ the replacement z —)• Cm in the expression of the function / produces 
constants due to the fact that / verifies the equation Polubarinova Galin, as will be 
proved further below. These constants are denoted /3m- After that we take the imaginary 
part 


'TV+l 

. 1=1 


x-Ci (t) 

1 


dCi (t) 


+Im 

N+l 

> of _ 

^ x-Ci (t) dt 

V+l N+l ^ 

aiak 

1=1 k=l 


dCk (t) 


x-Ci{t)x- Ck (t) dt 


= 1 


B. The constraints resnlting from the P-G equation 


1. The algebraic system 


(A.12) 


The invariants of the solutions to the equations P-G are 

N+l 


fdk = fk (Cfc, t) = Ck-it-i^ai log {Ck - Cl) 

Here we apply the operator d/dt 


1=1 


N+l 


„ dCk (t) . . 

0 = - t :- Z — Z 2^ 011 


dCk (t) dCi (t) 


dt 


Ck {t) - Cl {t) V 


dt 


(A.13) 


(A.14) 
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Or, write first the real and imaginary parts of the constants I3k- We have 

log (Cfc - 0) = log [(^fc -^i) + i {-Vk - Vi)] 

= log(|(^fc -^z) -i{Vk + Vi)\) 

+i arg [(^fc - 6) - * iVk + Vi)] 


The modulus is 


log[|(^fc-6) -i{Vk + Vi)\] 

= \ log [i^k - + iVk + Vif] 


The phase is 

arg [(^fc -^i) -i iVk + Vi)] = arctan 

When the abscissa — 6 is negative the complex argument Cfc — Ci crosses the cut 
have to add in. 

In addition, when two singularities are such that Cfc = Cz we have 

Cfc - Cz = 0 Cfc - Cz = 0 and r]k + r]i-)■ 2T]k 

. , [Vk + m 

we have arctan I --— 

V Cfc - Cz 



/ Vk + Vi \ 

Uk-^iJ 


Then 


l^'k = ik-^a[ arctan 


l^k 


■-Trai. — TT 


Vk + Vi 

Cfc Cz 


a[ + al log {2irik 


Zgfe —?!<0 


+x Y [(Cfc - iif + (hfc + vif] 


l^k 


and for the imaginary part 


l^k = -t-Vk- a'fc log {2r]k) - - tt Y “z' 


Zgfc—6<o 


— 'Y^ ^'i arctan 


l^k 


Vk + Vi 

Cfc Cz 


-^Y^'i [(Cfc - Cz)^ + ivk + vif_ 


l^k 


(A.15) 


(A.16) 


(A.17) 

and we 


(A.18) 


(A,19) 


(A.20) 
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Note that the in-determination induced by the function log is made explicit by the mul¬ 
tivalued function arctan. The arguments of the arctan are real and the hrst determination 
Cretan is hnite for any choice of variables. 


The unknown functions are the derivatives with respect to time, of the real part and 
imaginary ( 77 ^) functions of time-dependent positions of the singularities. 

dik dT]k 


dt ’ dt 

Let us introduce systematic standard notations. We note 
A- (1) ^ A (2) ^ A (A + 1) ^ 

A (A + 2) = ^, A (A + 3) = A (A + 1 + A + 1) = 


The variables 


X (i), for i = 1, 2N -|- 2 


(A. 21 ) 


(A. 22 ) 


(A.23) 


are the unknown. 


2. The coefficients as result from the equation for the real part (d/3ydt = 0) 


We prepare the coefficients in the linear system for ^ and ^ by calculating the deriva¬ 
tives 


d ^ fVk + m 

— arctan --— 

dt \^k - 6 


(A.24) 


l+(i3 




fd^ 

1 

Si- 

' \ dt 

dt J \ 


(6 - 6 )^ + iVk + Vif 


-ii\dt dt J - ^lY \ dt dt 


u c \ i^dk ^ drii\ .[dik dii 


Further 


^ log [(6 - iif + (% + mf] 


(A.25) 


{ik — ^l) + iVk + Vl) 


Ofc CM ,07 I ^ fdVk dqi 

+2fe + ,„) —+ - 


and 


^ log {2r]k) = ^ [log {r]k) + log ( 2 )] 
1 d7]k 


(A.26) 


Tjk dt 
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(A.27) 


Then we write 

dP'k 

dt 

_ d^ 
dt 

_ I _ ^ _ 

^ + iVk + Vif 

T]k dt 




d^k _ d^ 

dt dt 


E 


Of; 


l^k ~ + (dk + Vl) 

In Eq. k. For k = 1, N + 1. 
Coefficient of ^ is 

1 + 




(A.28) 


l^k 

l^k 


- - ..2^. ^ ^2 [- iVk + Vl)] 

(Cfc - 6) + iVk + Vl) 


[(6 - 6)] 


(Cfc ~ Ci) + iVk + Vl) 


Coefficient of ^ , for / = 1, A^ + 1 but I ^ k; it is 


~cy.] 


[iVk + Vl)] 


^ (6 - 6 ) + ivk + viY 

JZ ^ x2 I / ~ ^ ~ ^l)] 

{^k - 6 ) + iVk + Vl) 


(A.29) 


Coefficient of this means the the subscript of rji is the same as the number of the line 


k. The coefficient is 


E“: 

l^k 

1 

^Vk 

E-" 

l^k 


I " 

+au — 


0^7 


(C C ^2 / n2 ^l)] 

(Cfc - ^l) + iVk + Vl) 


[iVk + Vl)] 


(A.30) 


(Cfc - 6) + ivk + viY 


Coefficient of for / = 1, A^ +1 with the constraint I ^ k; the subscript of rji is different 
from the number of the line k. The coefficient is 

1 


-a 


[(6 - 6)] 


(Cfc - 6) + ivk + viY 

+0^1 Tc — ^ s2 [ivk + Vl)] 
{^k - 6) + iVk + Vl) 


(A.31) 
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Free term; it is 0. 


a. The list of types of coefficients. We note that there are few types. 




for ^ ,l = l,...,k-l 


for ^ ,l=k + l,...,N + l 
for rj ,1 = l....,k-l 


for 77,1 =k + l,...,N + ] 
0 for 1 =k 

O for , 1 = k 


Figure 15: The schematic representation of the matrix of the system. 


First divide the full {2N + 2) x {2N + 2) matrix into four squares each of (A^ + 1) x 
(N + l). 

There are 6 regions along a line of the matrix A{I, J). 
the number of types is however less due to repetitions. 

We note that the coefficients for line I , coming from ^ but for J 7 ^ J is one type. They 
are to be found for “columns” J that are in the first Jacobi square matrix 

l<J<iV + l, l<J<iV + l (A.32) 


The first part is 


l<J<iV + l, 1< J < I 
This is the type W1 {k, 1). The expression is 


(A.33) 


Wl{k,l) = 


(A.34) 


-a 


-{-Q. 


+ ivk + viY 

, _ 1 _ 

^ ifk - + iVk + Vif 


[ivk + Vi)] 

[- (Yk - 6 )] 
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where 


k = I 
I = J 


(A.35) 


Then we have a different expression for the local diagonal J = I. 


J = I 


This type is hT2 {k, k). The expression is 
W2{k,k) = 


1 + 


l^k 

E 

l^k 


Of; 


7-- ^.2^. ^ .2 [- iVk + Vl)] 

{^k - ^l) + iVk + Vl) 


where 


(6 - 6 ) + iVk + ViY 


k = I = J 


(A.36) 


(A.37) 


(A.38) 


Along the line, for columns less than J = N + 1 which is the end of the hrst square, 


1</<A^ + 1, J + 1<J<A^ + 1 


the type is again W1 {k, /). 


W1 (k,l) = 


-a 


Ter 


where 


^Yk - Yif + (Vk + ViY 
1 

{Yk — Yl) + iVk + Vl) 


k = I 
I = J 


[iVk + Vl)] 

[- (Yk - ^z)] 


(A.39) 


(A.40) 


(A.41) 


Along the line J, we now go to the second square 


1<I <N + 1 , (ArTl)Tl<J<(ArTl)TArTl 


(A.42) 
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These columns come from Here again we have to make the difference between J ^ I 
and J = I. 

For the first columns, before the local diagonal, 


l</<iV + l, (iV + l) + l<J<(iV + l) + / 


(A.43) 


the type is hF3 (/, J). The expression is adapted to this range of indices. 


W3{k,l) = 


—a 


^ (Cfc - + iVk + Vif 


[(6 - 6 )] 


+«"77- ^ s2 [iVk + Vl)] 

{^k - 6) + iVk + Vl) 


with 


(A.44) 


k = I (A.45) 

/ = J - (iV + 1) 


Then comes the local diagonal in this square matrix 


1 < / < A^ + 1 , with J - (A^ + 1) = / 
for {N + 1) + 1<J<{N + 1) + N + 1 


(A.46) 


This is type hF4 [k, /). The expression is 


where 


WA {k, k) = 


l^k 

n 

^k 


(C C\^ A I I ^2 '^0] 

- il) + {Vk + Vl) 


A-al — 


Vk 


l^k 


(Cfc ~ ^l) + iVk + Vl) 


[iVk + Vl)] 


k = I = J-{N + 1) 


(A.47) 


(A.48) 


Next along the same line I, 


1<I <N + 1 , (Ar + l) + / + l<J<(Ar + l) + Ar + l 


(A.49) 
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the type is again W3 {k,l). The expression is 


where 


W3{k,l) = 


—a 


^ i^k - + iVk + Vif 


[(6 - 6 )] 


+«"77- ^ s2 [iVk + Vl)] 

{^k - ^l) + iVk + Vl) 


k = I 

l = j-(N + l) 


for the range (J, J) written above. 


(A.50) 


(A.51) 


Let us consider application to few lines. 
b. The line for the unknown X (1). This is 

The first case consists of: the indice (/) of the coefficient coincides with the number of 
the line (which is A; = 1). 

This hrst case is the diagonal element. It comes from The model is W2{k,k) for 
A; = / = 1. 


A (1,1) = 


1 + 
N+l 

1^1 

N+l 


E 

1^1 


Ca.] 


- - ^.2^. ^ s2 [- (ffi + Vl)] 

(6 - 6) + (hi + Vl) 


[(6 - 6)] 


(6 - ^i) + (hi + viY 


(A.52) 


The other cases of entries coming from ^ have the indice of the matrix entry (coefficient, 
/) different of the number of the line (A; = 1). 

They are non-diagonal elements, upper to local diagonal (which in this particular case is 
the main diagonal); first N non-diagonal - upper - elements. They come from for 1^1. 
The model is hhl (J, J), for J = 1 and / -|-1<J<A^-|-1. 


A(1,J) = 


“"(6-0)^ + (hi + h.)^ 

- ij)^ + iVi + Vj)"^ 


[(hi + hj)] 

[- (^1 - O )] 


(A.53) 
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Next + 1 entries on the first line originate from 

The real position in the line of must be found by adding to its indice I of ^ the number 
of coefficients already considered, iV + 1. 

They are classihed according to the same criterium. The indice (/) is equal or not with 
the number of the line {k = 1). 

The first case is I = k = 1, the indice is equal to the number of the line. It comes from 
The model is WA {k,l), adapted for the second indice to be returned to the / = I. 


A{1,N + 1 + 1) = 

Oil — 

hi 


(A.54) 




E 


( 6 - 6 ) +{vi+viy 


[(6 - 6 )] 


// 

ai 


77- .^2^. ^ ^2 [(hi + Vl)] 

(6-6) +(hi + h6 


which means A{1,N + 2). 


The next N come from ^ . The indice / in NOT equal with the number of the line 
k = 1. The type is again W3 {k, 1) adapted for 

J = 1 and (iV + 1) + 1 < J < (iV + 1) + iV + 1 (A.55) 


with 


A (1, A^ + 1 + /) — 
1 


(A.56) 


—a'l 


+ 0 ;/ 


(c cA\( I A 
(6-6) +(hi + h6 


y'j-c - ^ [(hi + Vl)] 

(6 - 6) + (hi + Vl) 

ioT I = 2, N + 1, which means A (1, A^ + 3 ),A (1, 2N + 2). 


c. The line for the unknown X (2). The first elements are non-diagonal elements; here 
are the lower part. 

For this case the indice I is NOT equal with the number of the line k = 2. 
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They come from ^ for I = 1,2 — 1. This means a unique element. The type is 
W1 {I, J) adapted for J = 2 and J = 1. 


A (2,/) = 


-a 


-\-(y 


^ (6 - 6 )^ + iv 2 + viY 
1 

^ (6 - Ylf + (h2 + Vlf 


[{V2 + Vi)] 

[- (6 - 6)] 


for I = 1,..., 1, a single value. 


(A.57) 


Next entry consists of the case where the indice I of the entry (coefficient) coincides with 
the number of the line k = 2. 

The diagonal element. It comes from The type is W2 (J, J) for J = I = 2. 


A (2,2) = 

1 + 

1^2 

1^2 


- - ..2^. ^ ^2 [- iV2 + Vl)] 

(6 - Yi) + {V2 + Vl) 


a 


^ (6 - Yif + iv 2 + vif 


[(6 - 6)] 


(A.58) 


Next cases are still coming from the set ^ but their indice I does not coincide with the 
number of the line /c = 2, it is greater than k = 2. The type is again W1 (J, J) for J > I = 2. 

A{2,1)= (A.59) 


Their number is I 


—a 




'^Tc - .^2^. ^ ^2 [(^2 + Vl)] 

(6 - il) + {V2 + Vl) 

2 1 -(6-6)] 


^ (6 - 6)^ + {V2 + Vl) 


2 + 1,..., A^ + 1 which is / = 3, A^ + 1. 


Next entries (coefficients) are upper diagonal elements and come from 
We must make difference between (1) / = 1, ...,N + 1 but I ^ 2 and (2) / = 2 , where 
the index I only refers to The indices for the full matrix are calculated adding N + 1, 
corresponding to all the positions occupied by 

The first case consists of the entries coming from ^ with I less than the number of the 
line k = 2. The type of these terms is W3 (J, J), for J — {N + !)</. This group here only 
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contains one element, since we are at the line 2. Then I = 1. This coefficient comes from 

drji 

dt 


A{2,N + 1 + 1) = 

-«1 -- 9 [(6 - 6)1 

{^2 - + {V2 + Vi? 

+«177 - ..2^. ^ .2 [iV2 + hi)] 

(6-6) +(h2 + hi) 


(A.60) 


It follows the entry coming from ^ with the indice I coinciding with the number of the 
line k = 2, i.e. I = k = 2. It comes from The type is IT4 (/, I). 


A{2,N + l + l) = A{2,N + 1 + 2) = 
1 


// 
O^O 


A^+1 

E“i 

+2 

N+1 

E“" 

+2 


V2 


IC C+ \ I I \2 [^'^2 6)1 

(6 - 6) + (h2 + Vl) 


77- .^2^. ^ ^2 [(^2 + Vl)] 

(6 — 6) + (h2 + Vl) 


(A.61) 


The next group contains all other coefficients coming from ^ whose indice / does not 
coincide with the number of the line {i.e. k = 2). This means for I = 3,...,N + 1. The 
indice will be produced by adding I = 3,A^ + 1 to A^ + 1 (already taken by ^). The type 
is again W3 (/, J) for / = 2 and (A^ + 1) + 1 < J. 


A (2, A^ + 1 + /) = (only for / = 3,A^ + 1) 

-«(--^ [(6 - 6)1 


(A.62) 


^ (6 - 6)^ + (h2 + ViY 


+Q;/ 


(6 - 6) + (h 2 + viY 


[{V2 + Vl)] 


This hlls all entries on the line k = 2. The two exercises (for k = 1,2) can be used as a 
check for the computer code. 
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3. The coefficients resulting from the equation for the imaginary part 

(dl3'l/dt = 0) 


The time derivative of the imaginary part, fd'l- 


dt 

. drjk ,1 drjk 

= —i-:- Oi 


(A.63) 


dt 


'rjk dt 


E 

l^k 




(6 - 6 )^ + ivk + Vif 


1 


l^k 


(6 - 6 ) + ivk + viY 




In Eq. k. For k = 1, N + 1. 
Coefficient of ^ is 


-E 

l^k 




l^k 


- - ..2^. ^ s2 [- iVk + Vl)] 

{^k - 6) + iVk + Vl) 


[(6 - 6)] 


(Cfc — ^l) + iVk + Vl) 


Coefficient of For / = 1, iV + 1, but I ^ k. 


(A.64) 


~Ol.] 


[{Vk + Vl)] 


{ik — ^l) + iVk + Vl) 

^ 1 

(T \2 , / ~ ^ ~ 

[ik - ii) + {vk + Vl) 


Coefficient of it is 


(A.65) 


— 1 — 0)h — 

Vk 

-E 

l^k 


{ik — ^l) + iVk + Vl) 


- 6)] 


l^k 


(Cfc - 6) + ivk + viY 


[iVk + Vl)] 


Coefficient of it is 


(A.66) 


“0^7 


[(6 - 6)] 


^ (Cfc - 6)^ + ivk + viY 

1 

^ s2 [iVk + Vl)] 
{^k - 6) + iVk + Vl) 


(A.67) 
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for / = 1, iV + 1, but I 7 ^ k. 

The free term; it is, after transfering it to the right side 1 which comes from all equations 
of P''. 

a. The list of types of coeffieients as result from the time derivation of the imaginary 
part (dl3'f/dt = 0) There are 6 regions and the types are also only 4. 

Consider the line I, which must be 


(7V +1) + 1 < / < (iV + 1) + iV +1 (A.68) 

The hrst columns J are from ^ which are under the local diagonal. The type is hll (/c, /), 

ni{k,l)= (A.69) 

1 

^2 [iVk + Vi)] 


-01 : 


—a 


(6 - ^i) + iVk + ViY 
1 

i^k — ^l) + iVk + Vl) 


- i^k - ■C/)] 


where the indice k of the variables is I adapted by snbstracting the (A^ + 1) lines of the first 
Jacobi sqnare. 


k = I-{N + l) 
I = J 


(A.70) 


The range of the indices is 


{N+ 1)+ 1 < I < {N+ 1) + N+ 1 (A.71) 

1 < J < / - (Ar + 1) 


The element that is on the local diagonal has 


I-{N + 1) = J (A.72) 

It corresponds to ^ on the line /, with the property J = I — {N + 1). The type is Tt2 (J, J) 
with the expression 


n2{k,k) = 


-E“" 

l^k 

-E“! 

l^k 


(r J. ^2 [- iVk + Vl)] 

(Cfc - ^l) + iVk + Vl) 

TZ . n2 / I “ ^l)] 

(Cfc - ^l) + iVk + Vl) 


(A.73) 
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with 


k = I-{N + l) 


(A.74) 


The range of the indices is 


(AT + 1) + 1 < J < (iV + 1) + iV + 1 
J = J-(iV + l) 


(A.75) 


For columns that are beyond the local diagonal but still in the third square matrix, they 
come from 

at 

The type is again hll (/, J) with the expression 


ni {k,l) 


^( 6 - 6 )^ + iVk + Vif 
,_ 1 _ 

+ iVk + Vif 


[ivk + Vi)] 
[- i^k - 6)] 


The range of indices is 


(AT + 1) + 1 < J < (iV + 1) + iV + 1 
I-{N+1)<J<N+1 


which are translated into 


(A.76) 


(A.77) 


k = I-{N + l) (A.78) 

I = J 


Now we continue along the line no. / into the fourth square. 

For an arbitrary line /, the first columns in the fourth square come from They are 
of the type hlS (J, J) with expression 




“0^7 


[(6 - 6 )] 


^ (Cfc - + ivk + viY 

1 

^ ^2 [iVk + Vi)] 

{^k - 6) + iVk + Vl) 


(A.79) 
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where 


k = I-{N + l) (A.80) 

/ = J - (iV + 1) 


with the range of the parameters 

(Ar+1) + 1 < J < (iV +1)+ iV +1 (A.81) 

(Ar+ 1) + 1 < J < / 


The element on the local diagonal comes from ^ and the type is hid (/, J). the expression 


IS 


with 


and the range 


hid {k, k) = 


l^k 


-I / ^ 

T — ttfc — 


.,1 

'Vk 


(6 - 6 ) + ivk + viY 


- 6)] 


l^k 


7~C - ^ s2 [iVk + Vl)] 

[Yk - Yl) + iVk + Vl) 


k = I-{N + l) = J-{N + l) 


(A.82) 


(A.83) 


(AT + 1) + 1 < J < (AT + 1) + AT + 1 
J = I 


(A.8d) 


Finally we have along the same line (J) the group of columns that are upper the local 
diagonal in the fourth square. 

They come from They are of the kl3 (J, J) type again. The expression is 


m{k,i) = 


“0^7 


[(6 - 6 )] 


^ (Yk - Yif + ivk + viY 

1 

^ ^2 [iVk + Vl)] 

[Yk - Yl) + iVk + Vl) 


(A.85) 


dd 



The indices are 

k = I-{N + l) (A.86) 

l = j-(N + l) 

and the range of indices is 

(Ar +1) + 1 < / < (iV +1)+ iV +1 (A.87) 

I<J<{N + 1) + N + 1 


b. The free term in the linear system of equations The free term is 


F{I) = 


0 for 


1 < / < A^ + 1 


+1 for (A^ + 1) + 1 < / < (A^ + 1) + A^ + 1 


4- The numerical implementation of the time evolution of the positions of 
the singularities 

Assume that we work with A^ + 1 singularities (in the circular case, there are N singu¬ 
larities, in the inhnite case there are A^ -|- 1 singularities). 

ik + irjk for k = l,...,N + 1 (A.89) 

The equations that we intend to solve are of general form 

= (A.90) 

for /c = 1,..., A^ -|- 1. Putting together the real and imaginary parts {^k,Vk) as independent 
unknown variables, the system in the form 

^ A(I,J)X(J) = F(I) (A.91) 

j=l,2N+2 

for / = 1, 2A^ + 2 

The expressions of A (J, J) depend of {^k,Vk)- Then we can see the sequence 
1. we start with a set of variables 

{fk,Vk) , A: = l,2iV + 2 (A.92) 
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2. calculate the matrix 

A{I,J) 

and the free term 


3. solve the linear system 


F{I) 




4. integrate in time using Runge-Kutta method the system of equations 


and hnd the new set Vk) at the time t + 6t. 

5. hnd the current interface r(f) 

6. iterate to point 2. 


(A.93) 

(A.94) 


(A.95) 


(A.96) 


Numerical studies consists of the inversion of the matrix equation 


AX = F 


(A.97) 


followed by time advancement using the Runge Kutta method (the d02pdf NAG routine). 


C. The time evolution of the mapping function / {z,t) 


Integrating Eqs.(A.96) we determine the positions of the singularities (k (t), k = 1, A^+ 
1, as function of time. Now we calculate 


f{z,t) = z — it (A.98) 

v+i 

-i ^ [a'l + ia”) 

i=i 

X {log(|x -ii + iy- iyi\) + zarg {x-^i + iy - iyi)] 


It has been showrp310,E3 that: if at f = 0 all imaginary parts of the singularities of the log 
function are positive 

r]i > 0 for V/ (A.99) 
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and all constants ai are real and positive 


a; G R (A.100) 

tt; > 0 for V/ 

then the imaginary parts of the zeros yi and of the singularities r]i remain positive for all 
time t > 0. This preserves the holomorphicity of /. 

The interface line in the mathematical plane {z) is the straight line which separates 
the two half-planes—cx) < x < oo, y = 0 and it is mapped through the complex function 
f {z,t) into a curved line T (t) in the real plane of physical variables X = Ile[f {z,t)], 
Y = Im [/ (z, t)]. We can identify the curve that corresponds to the complex line z = x + i0 
by taking a; G R and y = 0 in the expressions of X, Y. 

The general formula for the Real part is 

Re[/(;^,t)] (A.lOl) 

= X 

Af+l 

+ Yl ^”2 ^ 

1=1 

N+l 

+ Y^'i {.x + iy-fi + iyi) 

1=1 

where we choose a'l = 0 for all I, then y = 0. 

Xo{x,t) (A.102) 

v+i 

= a; - ^ a; arg (x - 4- iyi) 

i=i 

Similarly, we have the general form 

Im[/(^,t)] (A.103) 

= y-t 

N+l 

1=1 

N+l 

+ Y {.x + iy-fi + iyi) 

1=1 
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and take a” = 0 and y = Q 


Yq (x, t) = -t 


N+l 


[(a^ - + {y- Vif] 


1=1 


(A.104) 


The interface is formally defined as 

r{t) = Yo{Xo,t) (A.105) 

which means that the variable x is eliminated between Xq [x, t) and Yq (x, t). 

The phase of the logarithm must be contained in 

0 < arg {x — + irji) < tt (A.106) 

because f]i > 0 for V/. There is nothing special with x traversing the point if we think in 
terms of the vector in complex plane based in the origin (0, 0) and pointing to {x — rji). 
When x — .^z = — |e|<0 the vector is almost aligned with the imaginary axis and the angle 
it makes with the abscissa (the phase of the argument of the logarithm) is slightly greater 
than 7r/2. When x moves to become greater than , a; — = |e| > 0, the vector is still 

almost vertical but the angle with the abscissa is slightly less than 7r/2. Therefore when x 
traverses the position the phase of the logarithm smoothly changes around 7r/2. When the 
argument of the complex number (a; — + irji) is calculated numerically, we must take into 

account that the function arctan = (tan)~^ is determined between —7r/2 and 7r/2. Consider 

? 7 i > 0 (which is allways true) (A. 107) 

X — = + |e| (very small) 

which means that x approaches from the right, such that x—^i > 0. Then arg (a; — + ir]i) < 

I which is obtained with the arctan function, since 

arctan (+CX)) = ^ (A.108) 

and we can write for this case 

arg (x — ^1 + irji) = arctan ( ——— ) (A.109) 


arctan 


Vi 


X 
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For the opposite case, where x approaches from the left, 


r/z > 0 

x-ii = -\e 


(A.llO) 


the value of arg (x — + irji) is again close to 7r/2, arg (x — + irji) > |, but in this case 

the function arctan gives 


arctan 



arctan (— 00 ) 


TT 

2 


(A.lll) 


and if we want to use arctan to obtain the arg then we must correct the result, by adding vr 


arg (x 


+ iVi) = arctan 


x-ii 


+ TT 


The two functions become 


(A.112) 


Xo (x, t) 

= X 


N+l 


a'l arctan 


i=i 




x-ii 


TT 




a. 


Yo (x, t) = -t 


N+l 


^ log [(x - ^z)^ + {y- Vif] 

l=l 


D. Note on the stagnation points 

We hnd that the calculation of the imaginary part 

N+l 

Yo = -t-Y^ a[-\og [(x - ^z)^ + {y- Vif] 

1=1 

makes Yq (x, t) strongly (linarly) dependent of time t. Then apparently the lines representing 
the interface are always substantially shifted one relative to the other and the stagnation 
pointM do not apper. However it can be proved that the linear term —t is cancelled by an 
opposite term arising from the sum in the second term. 
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The only possibility is that 




y = 0 this is the interface 

t 

rji -)■ exp 


a. 


This would be realized if the set of equations of time evolution of the would give 

- in the second half, for 77 ;, the general form 


A (/, J) 


dt 

dyi 

dt 


0 


or 


drji 

dt 


= - 7di + ••• 


a 


The diagonal of the matrix is 


hid (/c, k) = 


l^k 


-1 / ^ 

T — ttfc — 


.,1 

'Vk 


(6 - 6 ) + ivk + viY 


- 6)] 


l^k 


and we can approximate 


^ ^2 [iVk + Vl)] 

[Yk - Yl) + iVk + Vl) 




T — a',— — ... 

lyi 


We note that 


r]i ~ very small 

,1 

Of/ 

Vl 


-)■ 


> 1 


and the first term —1 can be neglected. Then 


dvi 1 

- 7 -- jVi + ••• 

dt a\ 
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and the solution is indeed 


m (t) 


( 0 ) 

rji exp 


a, 


(Here we see why must be > 0). 

We know that the singularities have evolved, qualitatively, in this way 


• the real parts are approaching one the other but this is a slow process. Approximately 
the initial positions are almost unchanged. 

• the imaginary part decreases rapidly, remains positive but approaches zero exponen¬ 
tially fast. 


Then we hnd that when we approach with x the position of one of the singularities 

a; - -)■ 0 


the term in the sum 

N+l 

^0 = -f - 5^ log [(x - + {y- riif] 

1=1 


becomes approximately 


Yn 


-t — a'l log 


viY 


where y 


log 

~ const > 0 


0 on the interface 


and this means that the tip [Xq it ), Yq (f)] for f —)■ cx), of the curve 

[6 (t) , r]i (f)] 

is almost hxed. These are the stagnation point^^. 

We note that the line Re [/ {z, t)] shows small but abrupt changes. These are not singu¬ 
larities however, but they are the result of the hnite precision of the representation of the 
two lines that are mapped by the conformal transformation. Essentially it is the degree of 
the detail when one approaches a singularity of the complex log function which introduces 
these jumps. 
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The shape of the Interface: Y=lmag(f) as function of X=Real(f) 



Figure 16: The time variation of the shape of the interface [X (t ), Y (f)] for the latest 11 
time moments, t = 90,100. The effect of imprecisions generated by {x — ^i, rji) —)■ (0, 0) 


consist of spurious lines that connect the hngers. They are made visible in Fig 20 


We explore the regions around the quasi-singularities with adapted mesh refinement. The 
spatial interval is x E [a^min, ^tmax] and we choose a number of average mesh intervals, M. 
The average mesh interval is 6x = (xmax — X min l/M. Now we choose a function / {x-,^i) to 
modulate 6x around a point x. We take 

NN 

f = c + ^aiexp 

1=1 

where c = constant, ai = amplitude of / for / = l,iViV = iV-|-l, the number of singularities 
Q. bi = half-width at inflection point. For example, for M = 1000, c = 1, a; = 40, 
bi = — ^i-i)/dienad di = constant factor = 10, the rehned mesh has 9700 intervals. The 

strongly non-uniform mesh allows a good precision in the regions where x — ^ ~ 0 but are 

not able to remove the imprecisions when also rji ~ 0. We conclude that special precautions 
must be taken for the late phase of the time evolution since in this regime rji becomes very 
small. In Figsj^ [^andj^we have avoided this region. 


- 6 )^ 

262 
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The lengths of the intervals on the refined mesh 



Figure 17: The lengths of the intervals of the rehned mesh on the line x = Re (z) in the 
“mathematical plane” . The non-uniform distribution is imposed by the need to explore 
carefully the regions around the quasi-singularities: (x — ^i) ~ 0, r/; ~ 0. 
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X=Real(f) as function of XL 



Figure 18: The X coordinate of the points on the interface {i.e. X = Re [/ (z)]) as 
function of a:, the coordinate of the abscissa on the mathematical plane z = x + iy. A 
small x—interval is plotted, to show the quasi-singular variation of Re [/ (z)] which is due 
to two situations where we have both x — 0 and rji ~ 0. These cases produce the 

errors of the interface profile. The result here is for t = 94. 
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Figure 19: 


Y=lmag(f) as function of XL 



5ame as Fig 16 but here Y = Im [/ (z)] is plotted. The quasi-singular behavior 

has the same origin. 
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The shape of the Interface: Y=lmag(f) as function of X=Real(f) 



Although a high degree of mesh refinement is used, the space variation of the interface is 
not correctly resolved around the two points shown in Figjl^and spurious lines are 
introduced, the same that are seen in FigjTdj 


Figure 20: A small region of the interface, from the same results as Figsj^and 19 
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B. APPENDIX. WRINKLED FRONTS AND CUSP SINGULARITIES 


For the examination of the cusp prohles, the equation of Sivashinsky type is solved in 
terms of a set of singularities (poles). The time dependence of the solution is encoded in 
the dynamics of the poles. The equations verihed by the poles ard^,h^ 

2N 


dz ■ 
dt 


cot 


Zk 


k=l,k^j 

+i^sign [Im(2;j)] 


(B.l) 


Here the poles are counted all, with the hrst N indices j = 1,..., A for poles and the last N 
indices j = N + 1,.., 2N for the conjugated poles. 


Zj+N — Zj 


(B.2) 


The equations are written for the real and imaginary parts of the poles 

Zj (t) = Xj (t) + iyj (t) 


and 



= z/ 


N 

E 

k=l^k^l 


Sin [Xi 


Xk) 


1 

cosh {yj -yk) - cos {xj - Xk) 
1 

cosh {yj + yk) - cos {xj - Xk) _ 


= \ _ sinh {yj - yk) _ 

dt k^^j ~ ~ ~ 

_ sinh {yj + yk) _' 

cosh {yj +yk) - cos {xj - Xk) _ 

+ 1 / coth {yj) 


-L 


These are the equations that are solved numerically. 


(B.3) 


(B.4) 


(B.5) 
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C. APPENDIX. A NOTE ON THE DISCRETE MODEL FOR THE 
BREAKING OF A RISING CONVECTIVE COLUMN 


The coupled lattice map model that is used to represent the physical process of phase 
competition is implemented numerically. Essentially the discrete nature of the problem (in 
particular the representation of the Laplacian operator) induces a certain stability of the 
phases, examined by Oppo and KapraP^. On a two-dimensional square lattice we initialize 
the held in one of the phases and add a small amplitude noise with smooth prohle. They 
are perturbations of Gaussian 2D shape, with positions, widths and amplitudes generated 
randomly. We then start the iteration and note the progress of the phase II into the region 
occupied initially almost completely by the phase I. 

One easily see formation of a spatial oscillation pattern whose characteristic extension 
is close to the unit cell of the lattice. This is connected with the choice of the diffusion 
coefficient and of the time advancement and are natural element of the model. We must 
however remove it if we want to use the evolving pattern of the phases to measure either 
the length of the interface or the area of the phases. This is a simple numerical operation 
but introduces a certain imprecision, which should not affect the application of this iterative 
discrete model to the problem of loss of compacity of rising convective columns. 

One can measure various quantities of interest, like the connectivity: how many compact 
patches of the initial convective column still exist in a horizontal plane; or, the area occupied 
by one of the phases and respectively, the length of the perimeter of the patches of one phase. 
For this we use the contour function (either Matlab or Fortran) and extract the set of closed 
curves that are defined by the same level. The dependence with the parameters of the 


iterative map can be studied, as shown in Fig. 14 
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